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^ • Abstract. Rediscovered by a systematic search, a forgotten class of integrable surfaces is shown 

to disprove the Finkel-Wu conjecture. The associated integrable nonlinear partial differential 
equation 



■ Zyy + il/z)xx +2 = 

(~| ■ possesses a zero curvature representation, a third-order symmetry, and a nonlocal transformation 

to the sine-Gordon equation (jy^^i — sin (j). We leave open the problem of finding a Backlund 
autotransformation and a recursion operator that would produce a local hierarchy. 
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Q>^ ' AMS classification scheme numbers: 53A05, 35Q53 

oo ■ 

O \ 1. Introduction 

(N 

^ . With this paper, we launch a project to classify integrable classes of surfaces. These 
are classes of surfaces whose Gauss-Mainardi-Codazzi equations are integrable in the 
sense of soliton theory. Our long-term goals include obtaining lists of integrable classes 
^ . as complete as computing resources permit, clarifying their mutual relations, and 
■ identifying known subcases. Our immediate goal is to demonstrate that the task is 
feasible and worth doing. 

The classical geometry of immersed surfaces in the Euclidean space is well known 
to be closely connected with the modern theory of integrable systems [35]. The Gauss- 
Weingarten equations of a moving frame always take the form 

= A^, = 5^. (1) 

where A, B are appropriate matrix functions. Integrability conditions of (1) are called the 
Gauss-Mainardi-Codazzi equations and take the form of a zero curvature representation 

Ay-B^, + [A,B] = {]. (2) 
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Equation (2) is invariant under a huge group of gauge transformations 

A' ^ S^S-^ + SAS-\ B' ^ SyS-^ + SBS-\ (3) 

induced by linear transformations ^' = 5'^' of the frame. Here 5" is an invertible 
functional matrix, which can be restricted to take values in the Lie group G associated 
with the Lie algebra Q matrices A, B belong to - typically 5o(3). 

The zero curvature representation (2) is the key ingredient in the soliton theory [15], 
where matrices A, B are additionally assumed to depend on what is called the spectral 
parameter. The essential requirement for solitonic integrability is that the spectral 
parameter cannot be removed by means of the gauge transformation (3). Consequently, 
if the matrices A, B can be modified so that they depend on a nonremovable parameter 
and still satisfy (2), then the corresponding Gauss-Mainardi-Codazzi equations are 
considered to be integrable in the sense of soliton theory, and their solutions are known 
as integrable or soliton surfaces [40] . 

Solitonic integrability can appear only when surfaces are subject to a constraint 
(such as being pseudospherical etc.). For numerous classical and recent examples see, 
e.g., the references [4, 35, 38] (or [16] in the projective setting). Workable tools to classify 
such constraints include all the general integrability criteria [31], which are, however, not 
immediately applicable to non-evolutionary systems [30]. Other methods take advantage 
of the already known non-parametric zero curvature representation (2), e.g., the method 
of extended symmetries by Cieslihski et al. [10, 11, 12]. 

In this paper we employ a recent method due to one of us [29]. Its essence 
can be summarized as follows: We attempt to extend the given non-parametric zero 
curvature representation (a seed) to a power series in terms of the spectral parameter. 
In the work [29], the relevant computable cohomological obstructions are identified. Two 
obstacles make this procedure not entirely algorithmic: The parameter-dependent zero 
curvature representation can exist only in an extension of the Lie algebra g and its jet 
order (the order of derivatives) can exceed that of the seed. If no obstructions are found, 
various ways exist to incorporate the true nonremovable parameter. 

2. Weingcirten surfaces 

To be of genuine interest in geometry, the determining constraint on integrable surfaces 
must be invariant with respect to coordinate changes. The general non-differential 
invariant constraint is a functional relation /(p, g) = between the principal curvatures 
p, q. Such a functional relation is characteristic of Weingarten surfaces, which have been 
a topic of continuous interest, especially in global differential geometry [20, 41, 38, 26] 
and computer graphics [8] . Well known to be integrable is the class of linear Weingarten 
surfaces [13, 35], characterized by a linear relation 

ak + hh + c = a,b,c = const (4) 

between the Gauss curvature k = pq and the mean curvature h = ^{p + q) (not to be 
mixed with a linear relation between the principal curvatures [23, 26]). Other integrable 
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classes of Weingarten surfaces that sporadically occur in the literature all have a 
differential defining relation (e.g., the Hazzidakis equation of the Bonnet surfaces [4, 5, 7]; 
a harmonicity condition of Schief's [37] generalized linear Weingarten surfaces) or the 
class is not determined by the functional relation f{p,q) = alone (e.g., [9]). 

So far, nothing contradicts the conjecture of Finkel [17, Conjecture 3.4] and Wu [43] 
that the only functional relation f{p, q) = to determine an integrable class of 
Weingarten surfaces is the hnear relation (4). Supporting arguments include Wu's [43] 
proof of non-existence of an so(3)-valued zero-curvature representation depending only 
on x-derivatives. Finkel's [17] argument roots in an unsuccessful search for higher- 
order symmetries and a (disputable, see [30, §2]) conjecture that integrability implies 
the existence of a local higher-order symmetry (actually the infinite hierarchy can be 
nonlocal, see also [31, §1.4.4.2]). 

Nevertheless, the main result of the present paper asserts that the simple relation 

= const (5) 

P Q 

between the main curvatures p, q, determines an integrable class of Weingarten surfaces. 
The associated nonlinear partial differential equation (21) has a parameter-dependent 
zero curvature representation (22) (outside the class considered in [43]), a third-order 
symmetry (24) (missed in [17]), and a recursion operator (25). 

Paradoxically enough, surfaces satisfying relation (5) were not unknown to nine- 
teenth century geometers. In view of their knowledge, our integrability result is not an 
entirely unexpected one. In fact, Ribaucour [34] established that the corresponding focal 
surfaces (evolutes) have a constant Gaussian curvature A; < (arc pseudospherical). 
Conversely, surfaces satisfying equation (5) arc involutes of pseudospherical surfaces. 
Moreover, the classical Bianchi transformation [2] is nothing but the induced 
correspondence between the two focal pseudospherical surfaces. Ribaucour's theorems 
are covered in Darboux [13] and early twentieth-century monographs, such as [3, 14, 18, 
42] . Later they became obsolete and forgotten as the induced Bianchi relation between 
pseudospherical surfaces became superseded by the classical Backlund transformation 
(the history is nicely reviewed by Prus and Sym in [32, Sect. 4]). 

The first examples of surfaces satisfying relation (5) also date to the nineteenth 
century. Lipschitz [25] derived a four-parametric family in terms of elliptic integrals. A 
particular subcase, the rotation surface of von Lilienthal [24], is the involute surface of 
the pseudosphere. 

The left-hand side of Equation (5) is equal to the difference of the principal radii 
of curvature at a point. This geometric quantity has a definite physical meaning, being 
associated with the interval of Sturm [39], also known as the astigmatic interval or 
the amplitude of astigmatism or simply the astigmatism [19] . A mirror or a refracting 
surface satisfying relation (5) will feature a constant amplitude of astigmatism in the 
normal directions. In the sequel, surfaces satisfying condition (5) will be called surfaces 
of constant astigmatism. Accordingly, the equation (21) to determine the surfaces of 
constant astigmatism will be called the constant astigmatism equation. 
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3. Preliminaries 

We shall consider surfaces r{x, y) parametrized by curvature lines. As is well known, the 
fundamental forms can be written as 

I = li^dx^ + 1;^ dy^, 

II = v?p dx^ + v^q dy^, 

where p, q are the principal curvatures. Coordinates y arc unique up to arbitrary 
changes x — X{x), y = Y [y). 

Let ^ — (ei, e2, n) denote the orthonormal frame, given by ei = Yx/u, e2 = rj/Zf, 
n = ei X e2. The Gauss-Weingarten equations 














( ° 




o\ 




V 






u 






% 







^ = 


_Vx 





vq 




V 








u 




v 


-up 


o) 






\ 


—vq 


0/ 



(6) 



are easily established. Their integrability conditions are the Gauss equation 



uUyy + vv^ 



V 



u 



-UtV 



UyVy + u^v^pq — 0, 



U V 

and the Mainardi-Codazzi equations 



(7) 



{p - q)uy + upy = 0, {q- p)vx + vqx = 0. 



(8) 



Consequently, the two so (3) matrices occurring in formulas (6) constitute a nonpara- 
metric zero curvature representation of the Gauss-Mainardi-Codazzi system (7), (8). 
Because of the isomorphism so(3, C) = sl(2, C), the same zero curvature representation 
can be alternatively written in terms of 2 x 2 matrices 



/ iUy 1 \ 

2v ^ 



\2 



1 1^2/ 
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(9) 



2w / 



Let us impose a constraint /(p, q) — 0. If nontrivial, it can be resolved with respect 
to one of the curvatures, say 



q = F{p), 



(10) 



which we assume henceforth. Then the Gauss-Mainardi-Codazzi system reduces 
substantially [8, 17, 43]. In particular, the Mainardi-Codazzi equations (8) have a general 
solution 



u 



-VqE', q = p 



E 



E' 



1 1 
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where E = E(p) is an arbitrary nonconstant function, E' = dE/dp, and and uq.vq are 
functions of x and y, respectively, removable by transformation x — J Uq dx, y — J Vq dy. 
Therefore, we can put uo — —vo — 1 without loss of generality, i.e.. 



1 , E 

«=^, v = E, q = p-—. 

The Gauss equation (7) then becomes 

Py, = E'E"p^^ + 2^pI + E\EE")'pI + EE'p" - E'p 



(12) 



Summarizing, the Gauss-Mainardi-Codazzi system of Weingarten surfaces reduces to 
the single equation (12). The classification problem considered in this paper is "for which 
choices of the function E{p) is the equation (12) integrable?" 

By substituting (11) into (9), we easily obtain a nonparametric zero curvature 
representation of equation (12), 



A, 



(]_Py_ 

2E' 
1 p 

\2E 



lp_\ 
2 E 

lPy_ 

'2E'J 



( i 



5n 



1 



EE"p^ 
2 2 

1 „v 1 



-{E'p-E) 
\-{E'p-E) --EE"p, 

which will be the starting point of the calculations to follow. 



(13) 



4. Cohomological criteria 

Readers not interested in details of the classification method can skip this section and 
continue to investigation of surfaces of constant astigmatism in Section 5. 

We use the formal theory of partial differential equations, which treats coordinates, 
unknown functions, and their derivatives as independent quantities. Equations can 
be conveniently represented as submanifolds in appropriate jet spaces [6]. All our 
considerations being local, we let J°° = J°°(M^, M) denote the space of oo-jets of smooth 
functions — )■ M. The base being equipped with coordinates the natural 
coordinates along fibres of correspond to p and its derivatives. These will be 

denoted p/, where / stands for a symmetric multiindex in x, y (including the "empty" 
multiindex such that p^ = p)- The usual total derivatives 

^ dx —'^^^ dpi^ ^ dy j'^^^ dpi 

can be viewed as acting on smooth functions defined on J°° (by definition, a smooth 
function locally depends on a finite number of coordinates). 

In J°°, we consider a submanifold Q determined by equation (12) and all its 
differential consequences obtained by taking successive total derivatives of both sides 
of (12). On Q, all derivatives of the form pjyy become expressible in terms of the others. 
Therefore, derivatives pi with y occurring no more than twice in / serve as natural 
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coordinates along the fibres of ^ — )■ R^. Being tangent to (12), tlie total derivatives 
admit a restriction to Q. We retain the same notation D^iDy for the restricted total 
derivatives. 

The essence of the adopted point of view can be summarized as follows: A function 
/ on satisfies f\g = if and only if / is zero as a consequence of equation (12). 
Prom now on we assume that all objects (hke the matrices A, B) are defined on Q. When 
writing 

{DyA-D,B + [A,B])\g^Q (14) 

we mean that the zero curvature condition (2) holds as a consequence of equation (12). 

In what follows, characteristic elements [27, 28, 36] play a crucial role. These 
are nonabelian analogues of characteristics of conservation laws [6]. For instance, the 
characteristic element of the initial zero curvature representation (13) is the s\{2, C)- 
matrix 



I'-- \ 

Co- . ^ . 

This immediately follows from the fact that 

DyAo - D^Bo + [Ao, Bo] = CqF, 
where 

F^pyy- E'E"p^^ - 2^pI - E^EE'JpI - p'EE' - pE\ 

so that the Gauss equation (12) can be written as F = 0. 

Let A = A{X), B = B{X) be the parametric zero curvature representation sought, 
C = C{X) the corresponding characteristic element. Besides (14), they will also satisfy 
the formula [27] 



0, (15) 



with I running over all symmetric multiindices, including the empty one. Here = 
Dj. — [A, ■], Dy = Dy — [B, ■], the other values being obtained by composition, which 
can be taken in any order since (14) implies that Dx,Dy commute. 

Characteristic elements of gauge equivalent zero curvature representations are 
conjugate (similar). This allows us to transform characteristic elements into the normal 
form with respect to conjugation, namely, the Jordan normal form. Since the matrix 
Co above is diagonal, it follows that for A sufficiently close to zero the characteristic 
element C(A) will be also diagonalizable. 



Integrability of Weingarten surfaces 



7 



However, diagonal matrices have a nontrivial stabilizer <S C SL(2,C) with respect 
to conjugation, which consists of diagonal matrices 



s 
1/s 



Gauge transformations from the group S (henceforth ^-transformations) preserve the 
characteristic elements C(A). Their gauge action on a general 5l(2)-valued zero curvature 
representation A, B is sufficiently simple: 




and similarly for B. Using iS-transformations, one can achieve a unique normal form 
of matrices A, B as follows: If ai2 7^ 0, then by setting s = (021/^12)^''*^ we turn 
A into a symmetric matrix, while in the remaining case ai2 = the zero curvature 
representation degenerates to a pair of conservation laws [28]. In other words, being 
symmetric is a normal form of nondegenerate zero curvature representations with respect 
to ^-transformations. 

Turning back to our original problem, we see that Bq is symmetric, and therefore 
the nearby matrices B(X) can also be symmetrized by an iS-transformation. A simple 
calculation shows that, by assuming diagonality of C(A) and symmetricity of B{X), we 
make the system (15) determined, hence solvable (actually, we fix the gauge). 

Summarizing, the computation of zero curvature representation has been reduced 
to solution of the determined system (14), (15) under a suitable choice of normal forms 
for C and B. However, this nonlinear system is still quite difficult to solve even with 
the help of computer algebra. To hnearize the system, the work [29] considers Taylor 
expansions 



^(A)-5^^,A^ 5(A) = ^5,A^ c{x) = j2CkX\ 



(16) 



k=0 



k=0 



k=0 



with Aq, Bq, Cq coming from the initial parameterless zero curvature representation (9). 
The condition of zero curvature for A{X), B{X) implies an infinite sequence of conditions 
of zero curvature for block triangular matrices 



/Ao 

Ai Ao 

\Am ' ' ' 



0\ 



Al Aq) 



51 



(Bq 

Bl Bq 

\Bm ■ ■ ■ 



■■■ 0\ 

■-. 

Bl Bq) 



(17) 



Characteristic elements C'™' assume the same form. Zero curvature representations 
^[m]^^H g^^g considered under the gauge group consisting of block triangular 



Integrability of Weingarten surfaces 



8 



matrices 



/ E 



0\ 



5'H — 



Si E 



m 



■•• 

Si e) 



with unit matrices E in the diagonal positions. By a cohomological argument presented 
in [29, Prop. 1], a nontrivial family A{X),B(X) with analytic dependence on A has 
expansions (16) such that Ai or Bi is not zero. 

Let (14) [™], (15) [™] denote the system obtained by substituting A A^"'\B -> 5H 
into system (14), (15), for arbitrary m > 0. Observe that systems (14)^^ (15)'™! are 
linear in their highest order unknowns ^4^, S^, Cm and can be solved sequentially. Then 
the applicable cohomological criterion can be summarized as follows. 

Proposition 1 ([29, Prop. 3]). Let m > 0. If Ai = Bi = for all solutions ^H^^H 
of system (14) t™', (15) ["^1, then there is no possibility to construct expansions (16) of 
order m, and consequently, the seed zero curvature representation Aq, Bq cannot belong 
to a nontrivial analytic family. 

Finally, to be able to solve system (14)["*], (15) '"^l, we need to know the normal 
forms of matrices A^'^\B^'^\ However, the normal forms for -B(A), C(A) established above 
immediately imply the same normal forms for (diagonal) and B^ (symmetric). 

5. Results 

In this section, we present the results of computation of the cohomological obstructions 
in the case of the nonparametric zero curvature representation (13) of equation (12). As 
a sub- result we obtain the first few coefficients A^, B^ of Taylor expansions (16). 

As we have seen in the preceding section (Proposition 1), the problem reduces to 
solving the system (14) '"^l, (15) '"^1 of linear differential equations in total derivatives, 
for increasing values of m. This is only possible under a suitable restriction on the jet 
order of the unknowns A^, Bj., Ck, k > 0. To start with, we assume dependence on the 
first-order jets at most. Upon expanding all total derivatives, equations (14)^™"', (15)^™"' 
become a large overdetermined system of linear partial differential equations. As such, 
the system is solvable by computing the passive (or involutive) form under a suitable 
(elimination) ranking [33]. 

Starting with m — 1, we checked that nonzero matrices Ai, Bi depending on second 
order derivatives exist for all possible determining relations (10). When incrementing 
m to 2, nontrivial conditions started to appear, but we also reached the boundaries of 
our available computing resources. Consequently, our present classification results are 
still incomplete. Nevertheless, we were able to obtain a passive system of differential 
equations in several cases. Moreover, in two cases we were able to find A2, B2 explicitly. 
One of them were the linear Weingarten surfaces (4). Their integrability is a well 
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established fact [35], the associated sine-Gordon equation (p^y — sin0 being a textbook 
example of integrability. The other class emerged as a solution 



E 



P 



c = const, 



of the ordinary differential equation 
E" _ /E'Y 2E' _ 1 _ 



Henceforth we concentrate on the solution (18). The coefficients u,v,q are easily 
found from (11) to be 



u 



p 



P + C 



pe 



l+c/p ' 



pc 
p + c 



The last equality shows that the condition of constant astigmatism (5) holds with the 
constant — 1/c on the right-hand side. The Gauss equation (12) becomes 



Pyy 



(? ^P + C n ^C^(C — p) n 



cp 



4 1+ 



P 



P 



2 1+- 



In principle, the cohomological method we applied can only prove nonintegrability 
and only indicate, but not prove, integrability. However, it was easy to guess an ansatz 
based on the form of and B^. By solving (14), (15) we obtained a A-dependent zero 
curvature representation 



A 



P P 



B 



ix + iy 

p 



^ 1+2- 

Ae p 

:^-v/A^e<'^9^ 
P 



\ 



P') 



-2 1 



'Px_ 
P' 



\/A^ + Ace"^ 



(19) 



V 



VFTa 



-1-2- 



ce 



-Ac 



P_y_ 



v/FTa 



c e 



which reduces to the initial Aq, Bq given by (13) when A = — |. The dependence on py 
explains why this class of Weingarten surfaces is missing in Wu's paper [43] . 
Upon substitution 



X 



P 



4c 



^ ~^ |^|l/4' y |^|3/4 

the Gauss equation (12) simplifies to 
^yy'^Qz) +2 = 0, 



21n2;-Fln|c| -4 



(20) 



(21) 
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1 + 2A z„ 
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(22) 



Let us remark that one can remove the x-derivatives from A and ^-derivatives from B by 
the gauge transformation (3), albeit at the cost of introducing an exponential dependence 
on the spectral parameter. In (19) and (22), the corresponding gauge matrix is 



-Ac/p g 







and 5" 



z'/' 
z-'/V 



respectively. For instance, the pair (22) becomes 



A' 



'\ ^J\^^\zy (A + 1)^-^ 



B' 



■ivWA% ^^W^\z-^-'^ 



~2 



Equation (21) has obvious translational symmetries d^, dy, the scaling symmetry 
2zdz — xdx + ydy, and a discrete symmetry 



X ^ y, y ^ X, z ^ 



(23) 



Computation reveals also two third-order symmetries of equation (21). One of them has 
the generator 



2 5 1 2 2 2 2 

-^Z {QZx ZZy^Zxx ~l~ (^"^o; ^ ^j^)(^2; ^^2/)-^a: 



K 
2 



r-2 -2: 



(24) 



y/'^xy 



■ Z Zx,^: 



x^xy 



+ -^Z^i^Zx ZZy)zl j^z{Zx ZZy)z^, 

where 

V (^^ - zZyY + 4/ . 
The other is obtained by conjugation with the discrete symmetry (23). 
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Moreover, A. Sergyeyev succeeded in finding a recursion operator for equation (21) 
in the usual pseudodifi^erential form 

-ZyD;' + z^D-^Dy + 2zD-^Dy (25) 

(unpubUshed). As far as wc could see, the operator generates only nonlocal symmetries. 
We leave as an open problem to find a recursion operator that would generate the 
third-order symmetry (24). 

Let us conclude this section with some easy geometric observations. First of all, we 
can put c — 1 without loss of generality. This can be always achieved by rescaling the 
ambient Euclidean metric and, if necessary, changing the orientation. 

Now, the symmetries of the constant astigmatism equation (21) have the following 
geometric interpretation. Translation symmetries are simply reparametrizations of the 
surface. The scaling symmetry (j)^: x — > c'^x, y — > e~'^y, z — > e^'^'^z takes a given surface 
r{x,y) to the parallel surface r{x,y) + en{x,y). This is not surprising since parallel 
surfaces obviously have equal astigmatism in the corresponding points. Finally, swapping 
the orientation is another symmetry, which can be identified with a composition of 
the discrete symmetry (23) and the rescaling 0i. Hence, the discrete symmetry (23) 
corresponds to the change of the orientation followed by taking the parallel surface at 
the unit distance. 



6. Relation to pseudospherical surfaces 

As already mentioned in the introduction, nineteenth century geometers knew of a 
simple relation between pseudospheric surfaces and surfaces of constant astigmatism, 
even though they did not find the latter important enough to be named. In this section 
we reproduce some of their findings and derive a nonlocal transformation between the 
constant astigmatism equation (21) and the famous sine-Gordon equation. Again, we 
put c = 1 for simplicity, meaning that the associated focal surfaces will be of Gaussian 
curvature —1. 

The forthcoming calculations are conveniently performed in terms of the variable 
z given by formula (20) or a new variable w related to z by 

^ = e^"'. (26) 

Then we have 

M = (w-l)e"'', v^^, p = -, q^-. (27) 

e w — 1 w 

and the discrete symmetry (23) becomes simply 

X ^ y, y ^ X, w ^ —w. (28) 

Given a surface C, recall that its evolutes (also known as focal surfaces) are the loci of 
the principal centres of curvature of C. Obviously, a generic surface £ has two evolutes. 
They interchange positions under the change of the orientation. 
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Proposition 2 (Ribaucour [34]). Evolutes of surfaces of constant astigmatism are 
pseudospherical surfaces. 

Proof. Let r{x, y) be a surface parametrized by curvature lines. We use the orthonormal 
frame (ei,e2,n), where 

ei=r^/u, e2 = ry/v, n = ei x 62. 

Then the two evolutes C and C" are given by 

, n ,, n 

r=rH — , r — r -\ — . 

p q' 

respectively. An easy calculation using the Gauss-Weingarten formulas (6) shows that 

/ Px f Py , 9^ 

p p \ p/ 



r. = --2n+(l--jr., 



= — 2^, 



the unit normals being 

n' = n" = ^. 

u V 

Now assume r{x, y) to be a surface of constant astigmatism. By applying the 
substitutions (27) we obtain the first fundamental form of the evolutes in terms of w: 

V = {wa; dx + Wy dyf + e"^"" dy"^ = dw^ + c"^"" dy^, 
1" = e^*" dx^ + {w^ dx + Wy dyf = e^^ dx^ + dw\ 

These are the well known pseudospherical metrics in terms of geodesic coordinates w, y 
and w,x on the first and the second sheet, respectively. □ 

For further reference we also compute the second fundamental forms 

II' = _e-^, dx' + || dy', II" = e'-wy dx' - |f dy'. 

Proposition 2 provides as with a couple of transformations from the constant 
astigmatism equation (21) to the sine-Gordon equation. To write them explicitly, we 
need to equip C and C" with the asymptotic coordinates ^,r), i.e., the fundamental 
forms have to be 

I' ^d^' + 2 cos (f)' d^ dri + drf, II' = 2 sin d^ drj, 

I" ^d^' + 2 cos 0" d^ drj + drf , II" = 2 sin 0" d^ drj. 

Here (f)' and (f)" are the angles between the coordinate hues on C and respectively. 
Using the previous expression of fundamental forms I', II' and I", II" in terms of the 
variable w, we easily see that ^, rj can be obtained by the "reciprocal transformation" [35] 

= I Viw^ + e'^^Wyf + e'^ dx + l Vie-'^'w^ + Wyf + e"'"' dy, 

(29) 

dr]^\ V(w^ - ^^'"wyf + e'" dx-\ ^/ie-'^^w^-Wyf + e-''" dy. 
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These formulas are valid on both sheets and reflect another property established by 
Ribaucour [34], namely that asymptotic lines on C and C" correspond. 
Then the angle (p' associated with the first sheet satisfies 



cos (I) — 



^/{w, + e'^w.r + e'"' V{w. - e'^wyf + e^- ' 



sm0' = ^ 



while the angle (f)" associated with the second sheet satisfies a similar set of equations 
related by the substitution (28). 

Proposition 3. Let z{x,y) be a solution of the constant astigmatism equation (21), 
let w — ^\rz. Determine function (f)' by formula (30), and new coordinates ^,r) by the 
reciprocal transformation (29). Then ^'(^,77) is a solution of the sine-Gordon equation 
(p^jj — sine/). 

Another solution of the sine-Gordon equation can be obtained by combination with 
the discrete symmetry (28). The other symmetries (translation and scaling) do not lead 
to essentially new solutions. 

Now, it is easy to check that the evolutes of surfaces of constant astigmatism are 
related by the classical Bianchi transformation. Indeed, the corresponding points r' and 
r" have a constant distance equal to 1/p — 1/q. The corresponding normals n' = r^/w 
and n" = Vy/v are orthogonal. Finally, being directed along the normal vector n, the 
line joining the points r' and r" is tangent to both surfaces C and C". These three 
properties characterize the classical Bianchi transformation. The Bianchi transformation 
is, however, superseded by the classical Backlund transformation [1], where the condition 
on the angle between the normals is relaxed from being right to being constant. 

7. Surfaces of constant astigmatism as involutes 

In principle, all surfaces of constant astigmatism can be obtained from solutions of 

the sine-Gordon equation as involute surfaces, see, e.g., Darboux [13, §802], Bianchi [3, 
§130-§150] or Weatherburn [42, Ch. 8]. Geodesic nets on pseudospheric surfaces fall into 
three classes: hyperbolic, parabolic, and elliptic [3, §102]. Of them only the parabolic 
geodesic nets lead to surfaces of constant astigmatism [3, §136]. 

Recall that the sine-Gordon = sin describes surfaces of constant curvature —1 
in the asymptotic coordinates ^, rj. By definition, 

I = d^'^ + 2 cos (fyd^drj + dr]'^, II = 2 sin d^ dr], 

which leads to the Gauss- Weingarten equations 

cos rt — r„ , . , cos r„ — re 
r^^ = ^— r^^ = sm0n, r^^ = — 0^, 

(31) 

cos <pr^ — Vjj cos r^ — r^ 
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Recall that coordinates X,Y ona. pseudospheric surface are called parabolic geodesic 
if the first fundamental form can be written as 

I = dX2 + e^^dy2. 

To find the transformation from asymptotic to parabolic geodesic coordinates, observe 
that d^'^ + 2 cos (j)d^dr] + drf = dX"^ + e^^ dY'^ is equivalent to the system 

X| + e'^F^ = 1, X^Xr, + e'^r^i;, = cos 0, X^ + e'^F^ = 1. 

This system can be rewritten as 

Xf — cos a, 5^ sin a, 

^ ^ _ (32) 

Xj^ = cos /3, Yjj = e ^ sin /3, 

and 

= a - ^. (33) 

In fact, (33) could be also (f) = /3 — a, which can be reversed by changing the orientation 
of the surface. The new unknowns a and P can be identified with the angles between 
the geodesies and the two asymptotic coordinate lines. 
The integrability conditions of system (32) are 

= — sin a, a^, — — sin /3, (34) 

or, in view of relation (33), 

^^ = -sin(0 + ^), ^^ = -0^-sin^. (35) 

These are already compatible by virtue of the sine-Gordon equation. Prom equations (32) 
we obtain 

sin B sin a ( cos /9 cos a \ ^ 

= — -. — r^f + -. — rrr,, ry = -. — ri"« + - — r^r, e . 



sin0 sin0 \sin0 sin0 

With respect to a given geodesic net, the involute surface f is composed of individual 
involute curves to the geodesies, based on one and the same orthogonal fine Y — const. 
Hence, 

f r + (a - X)xx, 

where a is an arbitrary constant. With the help of equations (31), the fundamental forms 
I, II of the involute surface r can be routinely computed in asymptotic coordinates. In 
particular, the unit normal is n = and 

i = (X^ - X + i)(l - cos 2a) di^ + (2X - l)(cos(a + /3) - cos(/)) did^ 

+ (X' - X + |)(1 - cos 2^) d^^, 
n = (X - |)(cos2q; - 1) + (cos(a + /3) - cos0) d^dr] 

+ (X - i)(cos2/3- 1)^77^ 
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Hence, the principal radii of curvature are X, X — 1. The Gauss-Mainardi-Codazzi 
equations of the involute surface hold as a consequence of the sine-Gordon equation, 
the two equations (32) on X and the system (35) on /3. 

To obtain the corresponding solution of the constant astigmatism equation (21), 
we have to reparametrize the involute surfaces by curvature lines. Let x, y denote the 
new coordinates. We choose x — Y and define y by the compatible system of equations 

= e^ sin a, = e^ sin /3. (36) 

A routine calculation shows that e"^^^^'^^ is a solution of the constant astigmatism 
equation (21). Summarizing, we have the following proposition. 

Proposition 4. Let 4>{^,r]) be a solution of the sine-Gordon equation = sin0. Let 
a, (3 be solutions of the compatible equations 

P^ — —sma, ar, — — sin (3, a — (3 — (f). 

Determine functions X, x, y by equations 

dX — cos ad^-\- cos /3 dr], 
dx — e"^ (sin ad^-\- sin /3 drj) , 
dy — e^(sin ad^-\- sin ^ drj) . 

Then the function e~'^^^^'y'> is a solution of the constant astigmatism equation (21). 

Example 1. Von Lilienthal's surfaces (involutes of the pseudosphere) . Published in 
1887, these surfaces seem to have fallen into oblivion. Recall that the pseudosphere is 
a surface obtained by rotating the tractrix around its asymptote. The meridians are 
geodesies of the parabolic type and therefore von Lilienthal's surface is obtained by 
rotating the involute of the tractrix (which itself is the involute of the catenary) . 

In geodesic coordinates X,Y, the "upper half" of the pseudosphere has a 
parametrization 



whose first fundamental form is dX"^ + e '^-^ dY"^ (differs by the sign of X from the 
canonical form used in the preceding section). Then 



parametrizes a rotational surface, for every real constant a. The surface is regular for 
all a < 0. Otherwise it has a cuspidal edge at X = a, which is a circle of radius e~". 



e'^cosF \ 
e"^ sin Y 



r = 



X > 0, 




,^ _ Vl - e-2^ I 
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Another singularity that occurs for every a > 1 is the intersection with the rotation axis 
at X — a — 1. Choosing the orientation so that the normal vector is 



n = 



/— e ^cos 
-e~^ sin Y 



-2X 



I 



I.e. 



n swaps orientation when crossing either of the singularities), then 



I = 



II 



e^^-l 



e^^-1 



and the principal radii of curvature are X — a and X — a-\-l. The corresponding solution 
of the constant astigmatism equation (21) is 



1 



z — 



X 



,2(a-l) • 




a = 0.7 a = l a = 1.2 a = 1.7 

Figure 1. A gallery of von Lilienthal surfaces 



Plane sections of von Lilienthal surfaces for various values of the parameter a can 
be seen on Figure 1. Besides the rotation axis, each picture shows the tractrix, which is 
the plane section of the pseudosphere, and its involute curve, which is the plane section 
of the von Lilienthal surface. 
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We finish this example with a short exploration of the behaviour at the limits of 
the definition domain. For X = oo the surface closes up at a point on the rotation axis 
at the height a — 1 + In 2, where both principal radii of curvature are infinite (the zero 
height is that of the cusp of the tractrix). For X — > the surface vertically approaches 
a horizontal circle of diameter |1 — a|. Two surfaces r{X, Y) and —i{X, Y) can be glued 
along this circle to form a single surface of constant astigmatism 1. For a — 1 both glued 
surfaces have a cusp here. 

8. Conclusions and discussion 

Among the still incomplete results of classification of integrablc Weingarten surfaces, 
we have identified a class originally introduced and investigated by nineteenth-century 
geometers. The class, which we propose to call surfaces of constant astigmatism, is 
governed by the equation 



For this equation we found an sl(2)-valued zero curvature representation depending 
on a parameter, a third-order symmetry, and a nonlocal transformation to the sine- 
Gordon equation = sin 0. We had to leave aside the problem of finding a Backlund 
transformation as well as a recursion operator producing a hierarchy of local symmetries. 

It should be stressed that the classification problem of integrable surfaces is far from 
being easy. An obvious reason lies in the abundance of integrability-preserving ways to 
derive one surface from another. Clearly, parallel surfaces, evolutes, and involutes of 
integrable surfaces are integrable. On the differential equation level, the corresponding 
notion is that of the covering [22]. The integrable classes of surfaces must be either 
closed with respect to taking derived surfaces or the derivation must map one integrable 
class into another. 
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